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Abstract. 

We consider the quantum correlations (entanglement and quantum discord) 
dynamics of two coupled spin qubits with Dzyaloshinsky-Moriya interaction influenced 
by a locally external magnetic field along z-direction and coupled to bath spin-i 
particles as independent non-Markovian environment. We find that with increasing 
D z and decreasing J z , the value of entanglement and quantum discord increase for 
both antiferromagnetic and ferromagnetic materials. Not that, this growth is more for 
the ferromagnetic materials. In addition, we perceive that entanglement and quantum 
discord decrease with increased temperature and increased coupling constants between 
reduced system and bath. But, strong quantum correlations within the spins of bath 
reduce decoherence effects. We discuss about type of the constituent material of the 
central spins that it can speedup the quantum information processing and as a result, 
we perceive that one can improve and control the quantum information processing 
with correct selection of the properties of the reduced system (J, J z , D z ). 



PACS numbers: 03.67.Hk, 75.10.Pq, 05.40.-a, 03.67.Mn 



Quantum correlations of Tuio-Qubit coupled to bath spin. 



2 



1. Introduction 

In the actual dynamics of any real open quantum system we have faced with the various 
type of interactions between system and its surrounding environment. These interac- 
tions which can lead to decoherence, will cause the transition system from pure quantum 
states to mixed ones and change the quantum properties, especially quantum correla- 
tions. Entanglement and quantum discord which are two different faces of the quantum 
correlations without classical counterpart, due to their notable features in developing 
the idea of quantum computers and other quantum information devices have attracted 
much attention in recent years. Since both of them have been realized to perform- 
ing quantum information tasks, the investigation of decoherence dynamics of them is 
an important emerging field [U |2j [3j SJ El [6] . Now, much attention has been paid to 
the quantum correlation in spin system, such as the Ising model [7] and all kinds of 
Heisenberg XX, XXZ, XYZ models [3 HOI Q3J. Both Ising and XXZ models can be 
supplemented with a magnetic term, the so called Dzyaloshinskii-Moriya (DM) interac- 
tion. Such antisymmetric exchange interaction arises from the spin-orbit coupling and 
has many important consequences. Recently, CHEN Yi-Xin and YIN Zhi [?] have per- 
formed an interesting investigation on thermal quantum correlations in an anisotropic 
XXZ model with DM interaction and shown that quantum discord is more robust than 
entanglement concurrence versus temperature T. They find that the characteristic of 
quantum discord is unusual in this system and this possibly offers a potential solution 
to enhance entanglement of a system. 

Just as other quantum systems, spin systems are inevitably influenced by their en- 
vironment, especially the spin environment. The coupling of spin systems with a spin 
bath often leads to strong non-Markovian behavior [] which has many physical im- 
portance in solid state quantum information processors, such as systems based on the 
nuclear spin of donors in semiconductors [T2| IT3] or on the electron spin in quantum 
dots [Hj. However the master equations describing the non-Markovian dynamics are 
rarely exactly solvable, but recently in Ref. [6] the authors present an exact calculation 
to obtain the reduced density matrix of two coupled spins in a quantum Heisenberg XY 
spin environment in the thermodynamic limit at a finite temperature. For spin systems, 
much attention has been devoted to the quantum correlations arising in spin chains 
at thermal equilibrium with their environment. In this sense, we consider the thermal 
quantum discord and entanglement within two spin qubit anisotropic XXZ Heisenberg 
model with antisymmetric DM interaction in the presence the bath spin-| particles as 
environment. Since the dynamics of the system qubit in the model we study is highly 
non-Markovian and hence it is not expect the traditional Markovian master equations 
commonly used, we follow the novel technique in Ref. [6]. In section 2 we introduce 
the model Hamiltonian describing two-qubit system with DM interaction coupled to an 
XY spin chain and derive an analytic formula for the exact solution of non-Markovian 
dynamics. In section 3 we review the concept of entanglement and quantum discord 
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and in section 4 we discuss the effects of the D z (the z-component of the DM inter- 
action), J z (degree of anisotropy), J (J > corresponds to the antiferromagnetic case 
and J < corresponds to the ferromagnetic case), T (temperate), g (the bath-system 
coupling constant) and g (the inner-bath-spin coupling constant) on entanglement and 
quantum discord, separately. Our results show that stronger DM interaction and week 
degree of anisotropy can decreases the disruptive effects of the environment. In addition, 
increasing the both temperature and the bath-system coupling constant will reduce the 
amount of the quantum correlation. Whereas the growth of the inner-bath-spin coupling 
constant or strong quantum correlations within the spins of bath reduce decoherence 
and improve the quantum correlations. By comparing the quantum discord with en- 
tanglement between two qubits in the model, we have seen that quantum discord is 
more robust than entanglement versus increasing temperature T and decreasing D z . 
Conclusions are then presented in Section 5. 

2. Model and solution 

The quantum system we consider consists of two spin-| anisotropic particles with DM 
interaction influenced by a locally external magnetic field along z-direction coupled 
to bath spin- 1 particles as environment. Here the DM interaction is a supplemented 
magnetic term arising from any interaction of a particle's spin with its motion, which can 
be represented by the form Dij.(Si x Sj). Where D is the DM vector coupling and 




the sum is over the pairs of spins. To see the effect of the DM interaction, we choose the 

— * 

^-component of the anisotropic parameter D. The corresponding Hamiltonian of total 
system H tot can be written as the sum of Hamiltonians for the system itself H s , the bath 
spin Hf, and the coupling between them H s i>. By considering interaction between the 
two anisotropic system-spin particles as Heisenberg XXZ model, with DM interaction 
parameter D z , the Hamiltonian of system is given by 



H s — e{S{ s + S% s ) + J Z S{ S S2 S + J{Si s S 2s + S ls S2 S ) + iD z (Si s S 2s — 5 , ls 5'^). (1) 



We restrict the interaction between the iV-components of the bath spin and between 
the system-bath to such that can be describe as Heisenberg XY model 



For definiteness, the s subscript refers to the system and b the bath spin. Where, the 
parameter e characterize the intensity of the magnetic field applied along the z-axis and 
D z the z-component of the DM interaction. The coefficients J z , J, go and g correspond 
to the real coupling constants with J z , J > for the antiferromagnetic case and J z , J < 
for the ferromagnetic case. By using the collective angular momentum operators 
J± = J2f=i Sfb an d the Holstein-Primakoff transformation as J + = a 1 (y/N — a) a) 
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and J_ = (yN — a^a) a with [a, a'] = 1, the Hamiltonians of the bath spin Hj, and 
interaction H s b can be rewritten as 

H b = goW(l - ^)a 




H sb = g{(S£ + S+)\l 1 - ^ a + (S u + S 2s )a^ 1 - ^}. (3) 

By considering the thermodynamic limit (i.e., N — > oo) at a finite temperature these 
above equations are reduced to 

Hb = 2g a*a, 

H sb = 9{(S+ + St)a + (S u + S 2 -)at}. (4) 

Observe that, the bath spin is reduced into a single- mode bosonic thermal field with 
non-Markovian effect on the dynamics of the our system. Since this thermal field will 
not remain in a thermal equilibrium state as usually assumed for an environment with 
very large degrees of freedom, therefor, the master equation approach can not be useful. 
Following the special operator technique introduced in [6], by tracing over the degrees 
of freedom of the environment from density matrix 

Ptot(t) = e- m °* t p tot (0)e m °* t 1 (5) 

we can catch the exact non-Markovian dynamics of reduced density matrix for the 
system at arbitrarily finite temperatures. We assume that the initial density matrix for 
the total system can be described by a pure and separable state as Pi O t(0) = p s (0) ® Pb- 
Where pb refers to the initial density operator of the single-mode bosonic thermal 
field which at thermal equilibrium is represented by the Boltzmann distribution as 
Pb = j?exp[—H b /kBT] with the partition function Z = (1 — exp[— 2g /kBT})~ 1 , where fee 
is Boltzmanns constant which we henceforth set equal to one. When two spin particles 
of the system is initially prepared in a normalized state with maximally quantum 
correlation as 

|^ s (0)) = «|00) + /3|ll), (6) 

it is easy to check that, the reduced density matrix for the system in the standard basis 
|00), |01), |10), |11) has the form as 

p s (t) = ^Ti h [\a\ 2 (e- iHtott \m)e- Hh/T {m\e iHtott ) + \f3\\e- iHtott \ll)e~ Hh/T {ll\e iHtott ) 
+ aP*(e- iHtott \m)e- HhlT {ll\e iHtott ) + a* P(e- iHtott \ll)e- Hb/T (00\e iHtott )]. (7) 

Calculation of the above equation could be more difficult due to mutual and internal 
coupled of the system and environment. One way of circumventing this problem 
introduced in [6] by converting the time evolution equation of the system under the 
action of the total Hamiltonian into a set of coupled noncommuting operator variable 
equations. Then by turning the coupled noncommuting operator variable equations into 
commuting ones via introducing a new set of transformations on the operator variables, 
the trace over the environmental degrees of freedom can be performed and the exact 
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reduced dynamics of the system can be obtained. 
For the case studied here, we can see that 

e- iHtott \00) = A{t)\00) + B(t)\01) + C(t)\10) + D(t)\U). (8) 

Note that, the coefficients A(t), B(t), C(t) and D(t) are functions of operators a and a) 
and do not commute with each other. The Schrodinger equation for e~ tHtott \00) is 

A( e -^|00)) = -iH tot (e- iH ^\00)), (9) 

which by replacing the Eq. (8), transform to 4 coupled first-order differential equations 
of noncommuting operator variables as 

j t A(t) = - i{{ J z - e + 2g <fa)A(t) + ga\B{t) + C(t))}, 

^B(t) = - i{(-J z + 2g a! i a)B(t) + ( J - 2iD z )C(t) + gaA(t) + ga*D(t)}, 
dt 

^-C(t) = - i{{-J z + 2g a jf a)C(t) + (J + 2iD z )B(t) + gaA(t) + ga*D(t)}, 
dt 

j f D(t) = -i{(J z + e + 2g Q a!a)D{t)+ga(B(t) + C(t))}. (10) 

From Eq. (8), the initial conditions are given by 

A(0) = 1, B(0) = C(0) = D(0) = 0. (11) 

Since, operator variables which appear in above coupled differential equations do not 
commute with each other, we cannot use the conventional methods of solving coupled 
differential equations. Fortunately, H tot is of an effective Jaynes-Cumming type and 
it can be block-diagonalized. Therefor, by finding the proper transformations of 
noncommuting operator variables, we can rewrite the Eq. (10) as the coupled differential 
equations of complex-number variables. In this model, we can reach our goal by using 
transformations as 

A(t) = e-^^-^A^t), 
B(t) = ae-* 90 ^- 1 ^ Bi(t), 
C(t) = ae-^^-^dit), 

D(t) = aae- 2 * 90 **- 1 * D^t). (12) 

Observe that, operator variables Ai(t), Bi(t), C\{t) and Di(t) are functions of n — a) a 
and commute with each other. Under these transformations the coupled differential 
equations of noncommuting operator variables change to the coupled differential 
equations of complex-number and commuting operator variables as 

j^it) = -i{(J z -e + 2g )A 1 (t)+gh(B 1 (t) + C 1 (t))}, 

j t B^t) = -i{-J z B 1 (t) + (J + 2iD z )C 1 (t) + gA 1 (t) + g(n-l)D 1 (t)}, 

jC^t) = -t{-J z C 1 (t) + (J-2tD z )B 1 (t)+gA 1 (t) + g(n-l)D 1 (t)}, 
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jD^t) = - i{(-J z + e - 2g Q )D 1 (t) + g^t) + d(t))}. (13) 

Notice that, if the total Hamiltonian cannot be block-diagonalized, for example, for the 
spin anisotropic particles with x- or y-component DM interaction, the operator method 
used here will then not apply to solve the problem exactly. By solving Eq. (13) in the 
usual way with initial condition ^i(O) = 1 and -Bi(O) = C*i(0) = Di(0) = 0, we can 
evaluate the time evolution for the initial two spins state of 1 00) . 

A similar analysis as above can be made if the two spins is initially prepared in |11) 
state. Let 

e -«W| n ) = A(i)|00> + B(t)\01) + C(t)\10) + D(t)\U). (14) 

For this case, the proper transformations of noncommuting operator variables have the 
form as 

A(t) = a^aU- 2t9o(h+l)t A^t), 
B(t) = o) e~ 2ig ^ +1)t B^t), 
C(t) = aU~ 2l ^ fl+1)t C 1 {t), 

D(t) = e- 2i9o(fl+1)t D 1 (t). (15) 
Insertion of these transformations in to Eq. (10) yield 

j t Mt) = -i{(J z -e + 2g )A 1 (t)+gh(B 1 (t) + C 1 (t))} : 

jB.it) = -i{-J z B 1 (t) + (J + 2iD z )C 1 {t) + g(n + 2)A 1 (t) + gD^t)}, 

jC^t) = -i{-J z C 1 (t) + (J-2iD z )B 1 (t)+g(n + 2)A 1 (t)+gD 1 (t)}, 

jD.it) = - i{{J z + e - 2g )D 1 (t) + g(n + l)(5i(0 + C^t))}. (16) 

By solving these coupled differential equations of complex-number and commuting 
operator variables in the usual way with initial condition Ai(0) = Bi(0) = Ci(0) = 
and -Di(O) = 1, we can evaluate the time evolution for |11). From the results of Eqs. 
(13) and (16), the reduced density matrix Eq. (7) in the representation spanned by the 
two-qubit product states |00) = |0)i <g> |0) 2 , |01) = |0> L ® |1) 2 , |10) = |l)i <8> |0) 2 and 
|11) = |l)i ® |1) 2 can be written as 

//On p u \ 

p 22 P23 



Ps(t) = 



P32 P33 
\p*U P uJ 



17) 



with 



1 oo 

Pii = y £( l«! 2 Mt) A[(t) + („ + !)(„ + 2)|/3| 2 Ax(t) 4(0 ) e- 2n90/T 



Z n=0 

Aigt oo 

PiA = L r T,<*rMt)DU- 2noo/T 

Z n=0 
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P22 = P23 = P32 = P33 = 

oo 

+ |cf n B x {t) B\(t) + (n + 1)|/3| 2 B 1 {t) B{{t) ) e~ 2n ^ T } 

n=l 

Pi , = \m 2 D 1 {t)D\{t){l + e-^ T ) 

oo 

+ J2( W? n(n - 1) D x {t) D\{t) + |/3| 2 D 1 (t) D\(t) ) e" 2 "^} (18) 

n=2 

Unfortunately, achieved analytic expression of solutions of Eqs. (13) and (16) has no 
compact form and, thus, we do not show it here explicitly. 

3. Quantum correlations 

The role of entanglement in developing the idea of quantum computers [15] and sending 
information in novel ways, such as quantum teleportation or quantum cryptography, 
has turned this intrinsic principle of quantum mechanics into one of the most prolific 
topics. To investigate the entanglement dynamics of the our bipartite system, we apply 
Wootters concurrence [16J. The concurrence can be calculated explicitly from the time 
dependent density matrix p s (t) of the two spins 

C(p s (t)) = max{0, Ai - A 2 - A 3 - A 4 }, (19) 

where the quantities Aj are the square roots of the eigenvalues of the matrix d = 
p s {t) {o~y ® o~ y ) p*{t) [a y <S> c y ) , arranged in decreasing order. Here p* s (t) means the complex 
conjugation of p s (t), and a y is the Pauli matrix. The concurrence varies from zero for 
a separable state to one for a maximally entangled state. The quantum state Eq. (17) 
are entangled if and only if either P22P33 < \Pu\ 2 or P11P44 < I/O23 1 2 - Both conditions 
cannot hold simultaneously [17J. The entanglement of this state is obtained as 

C(p s (t)) = y/pnpu + \pu\ - \y/puPM - \pu\\ ~ 2p 2 2- (20) 

However, concurrence is not the only type of quantum correlations. The different kind 
of quantum correlations than entanglement, so called quantum discord, has interesting 
and significant applications in quantum information processing. Quantum discord is 
not always larger than entanglement [T8l [T9] . This indicates that discord is not simply 
a sum of entanglement and some other nonclassical correlation. The relation between 
quantum discord, entanglement, and classical correlation even for the simplest case of 
two entangled qubits is not yet clear. In a bipartite quantum state with density matrix 
operator p s (t), that is inclusive of two part and has the form as two-qubit X states [20J, 
quantum discord is introduced as the difference between the total correlation and the 
classical correlation with the following expression 

D(p s (t)) = I(p s (t))-C(p s (t)). (21) 

There 

I{ Ps (t)) = S( Pu ) + S{p 2s ) -E7;log 2 7;, (22) 
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is the quantum mutual information. In this notation 7^ are the eigenvalues of the density 
matrix p s (t) 

71 = ^(Pn + P44 + y/(pu -P44) 2 + 4|pi 4 | 2 ), 

72 = ^(pn + P44 - y/(fhi -P44) 2 + 4|pi 4 | 2 ), 

73 = ^(P22 + P33 + \[{P22 - P^Y + Mp23 1 2 ) , 

74 = ^(P22 + P33 - \/(P22 -P33) 2 +4|p 2 3| 2 ), (23) 

and 

SW) = - [{pll + p22) log 2 (pil + P22) + (P33 + P44) log 2 (p 2 2 + P33)], 

S{p2s) = ~ [(Pll + P33) log 2 (pil + P33) + (P22 + P44) l0g 2 (p 2 2 + P44)], (24) 

are the von Neumann entropy of p ls and p 2s , the marginal states of p s (t). In order 
to computing the classical correlation C(p s (t)), we must perform a suitable projection 
measure on subsystem 2. After this, the state p s {t) will change to the ensemble {pi',Pi} 
(i = 0, 1) with 

Pi = -(I®V\i){i\V*)p s (t)(I®V\i){i\V*), 

Pi 

Pt = Tr[(I®V\i){f\V*)p a (t)(I®V\i){i\V*)]. (25) 

Where V G SU(2) is a unitary operator with unit determinant. The ensemble {pi',Pi} 
can be characterized by their eigenvalues as 

«±(Po) = \ (1 ± —\/[(Pll ~ P33)k + (p 22 - P44^] 2 + © ), 

2 Po 

1 1 

P ± ( Pl ) = - (1 ± -v/[(pii-p 3 3)^ + (P22-P44)A;] 2 + e ), (26) 
z pi 

with 

6 = 4/c/(|p 14 | 2 + |p 23 | 2 + 2Re{p u p 23 )) - 16mRe(p u p 23 )) + 16n/m(p 14 p 23 )), (27) 
and their corresponding probabilities as 

PO = (Pll + P33)k + (p 22 + P44)/ , 

Pi = (Pn + P3s)l + (P22 + Pu)k . (28) 

The parameters m, n, k and / are dependent to the projection of the von Neumann 
measurement V on the Bloch sphere. By rewriting this operator as V = tl + iY.a with 
t, 2/1,2/2, 2/3 e # and t 2 + 2/ 2 + y\ + 2/3 = 1 we have 

m = (tyi + 2/22/3) 2 , n = (ty 2 - ym)(tyi + 2/22/3), k = t 2 + y\, I = 1 - k. (29) 

The classical correlation is obtained as C(p s (t)) = S(pi. s ) — min{y|(po S(po) +p± S(pi)). 
Therefore, to calculate the classical correlation and consequently quantum discord, we 
have to minimize the quantum conditional entropy (poS(po) + pi S(p\)) with respect 
to the von Neumann measurements. In our model the minimization above expression 
attained in k — I — 1/2 and m = n = 0. 
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4. Quantum correlations of two coupled spins in spin bath 

In this section, we will discuss the time variations of the quantum correlations for an- 
tiferromagnetic (J > 0) and ferromagnetic (J < 0) in our model. Since the explicit 
expressions of entanglement and the quantum discord are very complicated, here we 
skip the details and give our results in terms of figures. 

In order to demonstrate the properties of different DM coupling parameter D z on quan- 
tum correlations in antiferromagnetic case, in Fig. 1(a) the thermal entanglement and 
in Fig. 1(b) the thermal quantum discord versus time are plotted for different values 
of D z . For both figures the coupling constants are J = 2.0 and J z = 1.0 and other 
parameters are e = 0.5, g = g = 1.0 and T = 2.0. From these figures, its obvious 
that the vanishing rate of the entanglement and the quantum discord decreases as the 
DM interaction parameter D z increases. This phenomenon indicates that, the decays 
of the quantum correlations due to interaction with spin bath can be compensated by 
tuning the DM interaction. In addition, it is easy to see that both entanglement and 
the quantum discord are recovered in the long time as we expected from non-Markovian 
dynamics. The comparison between these two figures shows that the quantum discord 
is more robust than entanglement under destroyed conditions. 

Similarly, Fig. 2 shows the characteristics of different D z on quantum correlations 
in ferromagnetic case (J = —2.0). Fig. 2(a) illustrates the thermal entanglement 
and Fig. 2(b) illustrates the thermal quantum discord versus time with J z = 1.0, 
e = 0.5, g = g = 1.0 and T = 2.0. In this condition, the quantum correlations decrease 
with gentle slope than antiferromagnetic case. As is indicated in Figs. 2(a) and 2(b) 
quantum correlations in ferromagnetic case is more robust than antiferromagnetic case 
versus increasing the values of the DM interaction. 

In Figs. 3 and 4 to exhibit the effects of the coupling constant J z on quantum cor- 
relations in antiferromagnetic (J = 2.0), the thermal entanglement in Figs. 3(a) and 
4(a), and the thermal quantum discord in Figs. 3(b) and 4(b) versus time for different 
values of D z are illustrated when e = 0.5, go = g = 1.0 and T = 2.0. In Fig. 3 the 
coupling constant J z has a lesser value than the same in Figs. 1 and 4. The compar- 
ison between these three figures shows that increasing the coupling parameter J z can 
enhance the fluctuations of the quantum correlations. It seems that increasing the J z 
will decrease the effects of the DM interaction. In one side, the curves of the higher 
value of J z are closer than of the lower ones. On the other side, the diminution of the 
both entanglement and the quantum discord in the long time scale for the higher values 
of D z is more in system with strong coupling. So as a result, the diminution of the 
quantum correlations can be compensated both by increasing the D z and by decreasing 
J z of the reduced system. 
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The effects of temperature T on quantum correlations are illustrated in Fig. 5. Where 
entanglement and the quantum discord versus time and temperature T with J = 2.0, 
J z = 1.0, e = 0.5, go = g = 1.0 and D z = 2.0 are depicted in Figs. 5(a) and 5(b), respec- 
tively. We found that in the very small temperature, both entanglement and quantum 
discord have their maximum value, but the increased temperature can decrease them. 
This decreasing of quantum correlations can be explained as follows. In the very small 
temperature, no excitation for the spins of the bath will exist and the bath is in a thor- 
oughly polarized state with all spins down. On the other word, the effects of decoherence 
induced by the bath are not strong enough to weaken the quantum correlations, but 
in the high temperature, the spin orientation of the bath becomes the very disorder, 
consequently effects of thermal fluctuation of the spin bath overcome on the effects of 
the quantum fluctuations (due to the D z interaction). With intensity decrease of the 
connection between bath and the reduced system, the role of mediation of the spin bath 
weaken, therefore two qubits of the reduced system can't exchange information with 
together and the value of quantum correlations reduced [21]. Not that, the quantum 
discord unlike entanglement never be zero. As a result, in the such as open quantum 
systems that the inner-bath-spin coupling constant is the strong, one can with right 
tuning J z and D z improve the both entanglement and quantum discord. 

5. Conclusions 

In summary, we have investigated the quantum discord and entanglement within two- 
qubit anisotropic XXZ Heisenberg model with antisymmetric DM interaction tunable 
parameter, such as D z , influenced by a locally external magnetic field along z-direction 
in the presence spin bath. We obtain that with decreasing J z and increasing D z , the 
value of entanglement and the quantum discord increase for the antiferromagnetic. In 
the same conditions for the ferromagnetic the quantum correlations decrease with gentle 
slope than antiferromagnetic case. In addition, we obtain that increasing the values of T 
(temperature) will reduce the amount of the both entanglement and quantum discord. 
Finally, we propose that in the such as open quantum systems that we can't reduce the 
devastating effects of T and g, we can with right tuning J z and D z or exchanging the 
constituent material of the central spins improve the both quantities of the entanglement 
and the quantum discord. 
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Figure 1. 



Figure Caption 

• Fig. 1(a): (Color online) The time evolution of the entanglement for initial two- 
qubit state \ip s (0)) = ^dOO) + |11)) for the antiferromagnetic case J = 2.0 and 
J z = 1.0. Other parameters are e = 0.5, g = g = 1.0 and T = 2.0. 
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Figure 2. 



Figure Caption 

• Fig. 1(b): (Color online) The time evolution of the quantum discord for initial 
two-qubit state \ip 3 (0)) = -^(|00) + |11)) for the antiferromagnetic case J = 2.0 
and J z = 1.0. Other parameters are e = 0.5, go = g = 1.0 and T = 2.0. 
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Figure 3. 



Figure Caption 

• Fig. 2 (a): (Color online) The time evolution of the entanglement for initial two-qubit 
state iV's(O)) = ( 1 00) + |11)) for the ferromagnetic case J = —2.0 and J z = 0.0. 
Other parameters are e = 0.5, go = g = 1.0 and T = 2.0. 
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Figure 4. 



Figure Caption 

FIG. 2(b): (Color online) The time evolution of the quantum discord for initial 
two-qubit state \ip s (0)) = ^dOO) + |11)) for the ferromagnetic case J = —2.0 and 
J z = 0.0. Other parameters are e = 0.5, go = g = 1.0 and T = 2.0. 
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Figure 5. 



Figure Caption 

Fig. 3(a): (Color online) The time evolution of the entanglement for initial two- 
qubit state \ip s (0)) = ^(|00) + |11)) for the antiferromagnetic case J = 2.0 and 
J z = 2.0. Other parameters are e = 0.5, go — g = 1.0 and T = 2.0. 
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Figure 6. 



Figure Caption 

• Fig. 3(b): (Color online) The time evolution of the quantum discord for initial 
two-qubit state \ip 3 (0)) = -^(|00) + |11)) for the antiferromagnetic case J = 2.0 
and J z = 2.0. Other parameters are e = 0.5, go = g = 1.0 and T = 2.0. 
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Figure 7. 



Figure Caption 

• Fig. 4(a): (Color online) The time evolution of the entanglement for initial two- 
qubit state ^(O)) = -^(|00)+|11)) for the ferromagnetic case J = 2.0 and J z = 1.0. 
Other parameters are e = 0.5, g$ = g = 1.0 and T = 2.0. 
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Figure 8. 
Figure Caption 

• Fig. 4(b): (Color online) The time evolution of the quantum discord for initial 
two-qubit state \ip s (0)) = -^(jOO) + |11)) for the ferromagnetic case J = 2.0 and 
J z = 1.0. Other parameters are e = 0.5, go = 9 = 1-0 and T = 2.0. 
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Figure 9. 
Figure Caption 

• FIG. 5(a):(Color online) The time evolution of the entanglement versus T for initial 
two-qubit state \i[) s (0)) = — 7^ (|00) + |11)) for the antiferromagnetic case J = 2.0 
and J z = 1.0. Other parameters are e = 0.5, go = 9 = 1-0 and D z = 2.0. 




Figure 10. 
Figure Caption 

• Fig. 5(b): (Color online) The time evolution of the quantum discord versus T for 
initial two-qubit state \ip s (0)) = ^(|00) + |11)) for the antiferromagnetic case 
J = 2.0 and J z = 1.0. Other parameters are e = 0.5, g = g = 1.0 and D z = 2.0. 



